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Abstract 

The amplitude of the roughness-induced energetic disorder at the metal/organic interface is 
calculated. It was found that for moderately rough electrodes, the correction to the electrostatic 
image potential at the charge location is small. For this reason, roughness-induced energetic 
disorder cannot noticeably affect charge carrier injection, contrary to the recent reports. 
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Effective injection of charge carriers is a major requirement for efficient and reliable perfor- 
mance of electronic organic devices (light emitting diodes, thin film transistors, and others).— 
Understanding charge injection is intimately connected to the structure of the metal/organic 
interface (or the interface of an organic material with another conductive material, such as 
indium tin oxide), and the detailed knowledge of this structure is important for other pro- 
cesses and applications. Recent experimental studies indicate that the Richardson-Schottky 
thermionic injection 



is a good starting point for the description of injection process in organic devices.— >2i^ Here j is 
the injected current density, (p is the height of the barrier at the interface, and E is an applied 
electric field. At the same time, measurements at low temperature show that the decrease 
of the injection current density is much smaller than the anticipated decrease according to 
the Richarson-Schottky model.^^ It was suggested that the reason for this discrepancy is 
the effect of energetic disorder in the organic material.— Usually, the calculation of the effect 
of disorder on injection is carried out using disorder parameters estimated from the charge 
transport data;^^ this means that these parameters describe the disorder in the bulk of 
the organic material. The experimental data clearly indicates, though, that in some cases, 
a surface dipolar layer is formed at the metal-organic interface.^^ It is very reasonably to 
assume that this layer has some degree of disorder and, hence, will provide an additional 
contribution to the energetic disorder at the interface.—"^ 

Recently, roughness at the metal/organic interface was suggested as the source of addi- 
tional energetic disorder, localized near the interface,- A calculation of the standard deviation 
of the disorder cr{zQ) for a point charge e located at distance zo from the mean plane of the 
weakly rough metal surface having profile h{x, y) may be carried out in the following way 
(we assume that the mean plane of the electrode is located at z = 0). Let us suppose that 
zq <^ I, where / is the surface correlation length. Then we can consider the electrode surface 
at the vicinity of a charge as a flat plane and treat surface deviation h{x,y) from the mean 
plane as a constant. The change of the image potential at the charge location due to the 
shift of the surface position by h is (in the first order in h) 




(1) 
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and <7{zo) is estimated as 



a^(zo) = e^(N(-o)0 = ^, (3) 



where the angular brackets denote an average over the ensemble of realizations of the surface 
roughness, h"^ = {h?{x, y)) is the roughness variance and e is a dielectric constant. The mean 
plane of the electrode is defined in such a way that {h{x,y)) = 0. Eq. Q is exactly the 
result of Ref. Isl, though obtained in a much simpler way. 

The important parameter for charge injection is the energetic disorder directly at the 
interface, i.e., in first several layers of organic transport molecules adjacent to the electrode. 
If for the very first layer zq = 6 A, ho = 3.5 A, and e = 1, then a{zo) = 0.7 eV (Ref. S). Yet 
the validity of Eq. Q for short distances is very dubious because of the basic assumption 
that organic molecules in any particular layer of organic material are situated at a constant 
distance zq from the mean plane of the electrode^ (see Fig. ^). The size of a typical transport 



molecule (8-10 A) 
20 - 50 nm, Ref. 



is small in comparison to the surface correlation length (typically, / 



10[ ). In this situation it is natural to expect that a better model of the 



interface is one where the molecules follow the electrode profile, and any particular layer is 
located at the constant distance zq to the actual surface of the electrode (Fig. ^p). We will see 
that in this model a{zQ) differs drastically from the estimation in Eq. (jH)). The calculation 
of cr[zQ) along the profile of the rough metal surface is the major difference between our 
paper and the paper of Rahman and Maradudin,— where the general expression for the 
mean image potential for a rough dielectric interface was obtained. 

The potential for a point charge located at the vicinity of a rough metal surface obeys 
the Poisson equation 

Aip = 6 {f - fo) (4) 

with the boundary condition 

We assume that the roughness is small h^/l <^ 1, and Eq. (jlj) can be treated via a pertur- 
bation theory approach. We are going to calculate the leading contribution only. A possible 
approach to perform this calculation is to transform to coordinates 2;new = z — h{x,y), so 
for new z the boundary condition is set for z = 0. In the new coordinates, the Poisson 
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FIG. 1: Two models of the arrangement of organic molecules (ellipsoids) at the rough surface of 
the electrode (solid line): a) first layer is located at constant distance to the mean electrode planej^ 
b) molecules in the first layer are located at constant distance to the surface of the electrode. 



equation takes the form 



(i)'+(il)'+i 



2 



(6) 



2 (at + Si) - fA^h = S if- fo) 



where A_l is a two-dimensional (2D) Laplacian and we assume zq = (0,0, ^o)- Note that 
in the new coordinates the condition zq = const is approximately equivalent to a constant 
distance to the profile of the electrode (with a small correction proportional to /ig/^^ and 
insignificant to our analysis). 

Let us try to find a formal solution as a series 

oo 

(p = ^Vn, ^n^OiK"), (^„|^^o = 0, (7) 
n=0 

Mr) 



e \r — zo\ e\r + zo\ 
The first-order correction is 



<^,(r) = j dnG{f,n)J{n), (8) 

here (^(r, ri) is the Green function for the Laplace operator with zero boundary condition 
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at z = 0, while the source term is 

J(r) = 2 f 1^1^ + l^f") + ^A^h 

^ ' \ axaz ox ayoz ay i az ^ 



e a 

e dzQ 



2 (ft + ft) + ^^-"1 . (") 
P{f) = + 

^ ^ \r-zo\ \r+zo\ 

Note that, in our case, the correction to ipo depends not on h{x, y) itself, but on its derivatives 
and vanishes for h{x,y) = const, as it should be. 
The Green function has the forn>i^ 

^(^"'^"i) = 4^/ dke-'^'^P-'''^Gkiz,z,), (10) 

where k and p = {x,y) are 2D vectors, and the Green function Gk{z,zi) obeys the equation 

^'^'^ - ec, = 6{z - z,), G,(0, z,) = 0. (11) 



dz"" 

The solution of Eq. (fTT|) is 

Gkiz,zi) = — — sinh exp(— /c2;+), (12) 

rh 

z+ = max(z, zi), z_ = min(2;, Zi). 

We are going to calculate the correction (jH} for r = zq only, because o"^(-2o) = {ip1{zo)). A 
simple but lengthy calculation gives for ipi 

(^o) = j dkidk2hj;^_j:^ exp [-{ki + ^2)^0] [cosh (ki - k2) ^0 - 1] , (13) 

and is a Fourier transform of h{p). The integral (|T3|) can be simplified further. Let us 
make a transition to new vector coordinates 'p={ki + A;2)/2, q = ki — k2- Then 

Vi{zo) = ^J^^^2 j dqg{qzo)h^, (14) 

where 

g(q) = J dp exp (-R) (cosh Q - 1) , (15) 
R=\p+^q\ + \p-lq\, Q = \p+^q\-\p-lq\■ 
Function g{q) can be easily estimated for g ^ 1 and g ^ 1. In the first case we can expand 
the hyperbolic cosine in the integral ()15|) in power series of q, thus obtaining 

g{q) ^ j dp (^^^ exp (-2p) = (16) 
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FIG. 2: Function g{q). 

In the opposite case {q 3> 1) analysis shows that the only significant (and equal) contribu- 
tions to Eq. (fT^ goes from p ^ ±| g, so setting p= | g + s we have 

g{q) ^ J dsexpi-2s) = |. (17) 

The general behavior of g{q) is shown in Fig. |21 

Finally, the variance of roughness-induced energetic disorder is 

a\zo) = {^lizo)) = J dqg\qz,)C,. (18) 

This equation is the major result of this paper. Here is the Fourier transform of the 
surface correlation function, which we define in a usual way (assuming spatially homogeneous 
roughness) 

{himPi)) = hlC{p- pi) (19) 
with C(0) = 1, so {h'^{p)) = /iq. For homogeneous roughness 

If zo ^ /, then we can replace g{qzo) by its limit value of 7r/2, and in this case 

.^(.0) « (21) 
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This result is equivalent to Eq. Q, but it is valid only far away from the rough electrode 
surface. The reason for the equivalence of Eq. (PTj) and Eq. Q is the need to cancel the 
leading term in Eq. (|14p in the old (physical) coordinate system. Indeed, as it follows from 
Eq. (HH) for Zq > / 

V'l (^o) = ^,+o(\\n[h], (22) 
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where VL is some integral operator. Image potential at the charge location is 
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(23) 
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(here z^^ denotes the distance to the mean plane of the electrode). This result means that 
in the old (physical) coordinate system the correction to the image potential in the first 
order in h decays faster than 1/ ^-Zo^'^j for large distances. This is not surprising because 
this kind of decay is possible only for h ^ const, which is not the case for where 
many uncorrelated domains of the rough surface contribute to the image potential. 

All these intricacies are not important for charge injection, where a relevant distance to 
the surface of the electrode is small. If Zq -C then 



pip 



dqq oc 



(24) 



2567r%2 J -^-^ " ^2/4 
The later estimation is valid if C(p) can be characterized by the scale / only, and the integral 
in Eq. converges for q ^ oo. If we assume a Gaussian correlation function 



£1 

2P 



which is a good approximation for indium tin oxide electrodes^, then for zq <^ I 



a'izo) 



(25) 



(26) 



Eq. ()2fij) is similar to Eq. (jH)) with the only crucial difference: Zo is replaced by /. For the 
roughest electrode, mentioned in Ref. with ho = 4 nm and / = 14 nm, we have at the 
interface a ^ 0.01 eV. 

Let us consider the case, when the integral ()24p does not converge for g ^ oo. This is 
the case of the fractal rough surface with a correlation function 



n -1 



;i + /t2/2)i+-' 



A = Aira 



1 - 



;i + k^.py 



(27) 



here < a < 1 (Ref. |13|). In fact, any fractal surface can be realized as an intermediate 
asymptotic only, for some spatial scale range, thus the proper cut-off kc is assumed in Eq. 
(I27|) . For a clear physical reason (discrete nature of a real metal surface) kc ^ 1/a where 
a is a typical interatom distance, while kj ^ 1. This means that Zq^c ^ 1 ioT the organic 
layers closest to the metal surface. Hence, we can still use the small-q asymptotic of g{q) 
and 



vi^ Jo ' 



^'i^^) Wr !t dq^^^^ ^ (28) 



256(2-a)7re2/4 • 

In the most favorable for large cr case a f« 

and for /iq = 5 A the value of cr becomes comparable with the bulk value of 0.1 eV only for 
fcc ^ 0.1 A~^. Such value for k^ seems to be unreasonably large. Scanning microscopy data 
indicate that typically the Gaussian correlation function (j25p is a good approximation for 
rough electrode surfaces at nanometer scale.— Fractal surfaces have been indeed observed 
in clusters formed by small metal particles, but the relevant spatial scale was very different: 



even the size of the individual metal particle was no less than 10 nm (Refs. Il4ljl5[ ): in this 
case kc < 0.01 1"^ 

In conclusion, we have found that the contribution from image forces to the roughness- 
induced energetic disorder is typically too weak to provide a noticeable effect on injection, 
contrary to the analysis in Ref. ^ 
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